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Abstract 

In this paper, a summary of the basic simulation 
parameters and results of a new study for the Geostationary 
Operational Environmental Satellite (GOES) is shown. The 
study for GOES involves the simulation of minor 
modifications to the current spacecraft, so that the relative 
performance of these modifications can be analyzed. The 
first modification studied requires the placement of a baseline 
inertial reference unit, such as the Dry Rotor Inertia 
Reference Unit (DRIRU-II) or the Hemispherical Resonator 
Gyro (HRG), onto the spacecraft. The imager/sounder 
assembly is currently used to obtain landmark and/or star 
observations in order to compensate for spacecraft motions 
and external disturbances through ground processing. The 
study utilizes the imager/sounder assembly as another attitude 
sensor for on-board attitude determination. Also, the addition 
of star trackers is used to provide precise attitude knowledge. 

Introduction 

The current (GOES I-M) spacecraft specification for the 
knowledge requirement is 112 jirad . This requirement is 

met through ground processing 99% of the time in the 
east/west direction and 95% of the time in the north/south 
direction. The spacecraft specification for the within-frame 
registration is 42 pxad . The current spacecraft uses an Earth 
Sensor Assembly (ESA) to provide roll and pitch information. 
Yaw knowledge is not sensed. However, yaw control is 
achieved through roll/yaw coupling. A set of gyros based on 
the Digital Integration Rate Assembly (DIRA) also is on the 
current spacecraft. However, the DIRA has an operational 
lifetime of 2000 hours. Therefore, the on-board gyros are not 
used for mission mode attitude determination and control. 

An outline of the remainder of this paper proceeds as 
follows. First, the simulation model for the gyro, the ESA, 
and the imager/sounder assembly are shown. This includes 
the simulation parameters used for Earth clouds and Earth 
radiance/gradients effects in the ESA, and non-repeatable 
errors in the imager/sounder assembly. Then, the simulated 
attitude sensor and gyro measurements are used in a Kalman 
filter for attitude determination. Results are presented for two 


cases: 1) using the ESA, and 2) using both the ESA and 
imager/sounder assembly. Next, results using a star tracker 
are shown. This includes simulation results with and without 
the addition of gyros. Finally, conclusions are stated based on 
the simulation results. 



In this section, a brief overview of the simulation 
parameters for the gyro model, the ESA model, and the 
imager/sounder model is shown. The true angular velocity is 
assumed to be modeled by [1] 

% 0) 

where CD is the true angular velocity, <5^ is the gyro- 
determined angular velocity, and £ is the gyro drift vector, 

b = r\ 2 ( 2 ) 

The 3x1 vectors, i] 1 and are assumed to be modeled by 
a Gaussian white-noise process with 

%(0} = Q i = l,2 (3) 

£ {n,('h^( r ')} = a 8 y8( J -0 ij = l2 (4) 

where 

fi = p 7 3x3 ° 3 * 3 ] (5) 

„ ^3x3 3x3 _ 

The DRIRU-II drift-rate noise and measurement noise 
characteristics are given by c u =2.15x10^ prad/sec 3 ^ 2 

and o v = 0.206 |irad/ sec 1 / 2 . The nominal motion of the 
spacecraft involves a rotation once per orbit about the 
spacecraft’s y-axis. Therefore, the nominal angular velocity 
is given by 
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0 


CD = 



( 6 ) 


where ©„ is the orbit rotation (727 x 10 5 rad / secj . 


The ESA measures the spacecraft’s roll and pitch angles. 
These angles are measured with respect to a moving Earth 
frame. The gyros provide attitudes with respect to an 
inertially fixed frame (e.g., GCI). Since the body rotation 
axis is about the spacecraft’s y-axis, the body measurement 
vector is given by [2] 




-sin(p)cos(r) 

sin(r) 

cos(/?)cos(r) 


(7) 


The imager/sounder assembly can measure stars in . 

23°E / W x21°N /S field of view, outside of the Earth limb. 
The orbit-attitude tracking system contains a catalog of bright 
stars visible by the imager/sounder which can sense three 
stars at 45 second intervals. For simulation purposes these 
stars are assumed to be found in different quadrants in the 
field of view. The imager/sounder star windows are 
staggered so that the data is acquired every 15 minutes. The 
imager/sounder measures the tangent of two angles, Pj and 
P 2 , resulting in a body vector given by [2] 




•^1+tan 2 P] +tan 2 p 2 ” 


tanp]" 
tanp 2 
1 


( 12 ) 


The imager/sounder “measurements” are obtained by 
using the following model 


where r and p are the scanner roll and pitch angles, 
respectively. The inertial reference vector is given by 

I e =A T (q)B e (8) 

where q is the true quaternion (obtained by kinematic 

propagation using the true angular velocity). The ESA 
‘"measurements” are obtained by using the following model 

p^p + v p +Wp (9) 

whore v p is a zero-mean Gaussian process with a 3o value of 
0.02 degrees, and w p represents the non-repeatable errors 

due to Earth cloud and Earth radiance/gradients effects. The 
non-repeatable error is assumed to be modeled by the 
following discrete process 

1/2 

w p (i + \) = Aw p (i)+L[l-A 2 ) g(i ) (10a) 


A = exp(-4 At B) (10b) 


where At is the sampling interval (0.25 seconds for the ESA), 
B is the bandwidth (for weather purposes, this set to about 
1/6 days), L is the lo amplitude (experience has shown that 
this is about 200 prad), and g is a zoo-mean normal 
Gaussian process. This same error model is applied to the 
Earth roll “measurement.” Since the roll and pitch 
measurements from the Earth sensor are small, the body 
measurements can be approximated by 


Ib 



( 11 ) 


tanPj =tanp, +v frl . +w bi , i = 1,2 (13) 

where v b is a zero-mean Gaussian process with a 3a value of 
28 prad . The non-repeatable error in the imager/sounder is 
assumed to be modeled by the following process 


X = 

' 0 

r 

x+ 

o' 


L-®n 

0 


i 


w b =[l 0]x (14b) 

where i] is a zero-mean Gaussian process. The standard 
deviation of T) is selected such that the output of w b has a 
3o value of about 200 prad . 

Simulation Results 

Fra: this part of the study, an investigation of the relative 
performance between using on-board gyros and without the 
use of gyros was examined. Fra the on-board gyro case, a 
standard Kalman filter with a gyro propagated model was 
used fra attitude determination. The simulations were run fra 
six cases, which include: 1) ESA only with no non-repeatable 
(NR) errors, 2) ESA only with NR errors, 3) ESA and 
imager/sounder (I/S) with no NR ESA errors and no NR I/S 
razors, 4) ESA and I/S with no NR ESA razors and with NR 
I/S razors, 5) ESA and I/S with NR ESA errors and no NR I/S 
errors, and 6) ESA and I/S with both NR ESA errors and NR 
I/S errors. 

The first two cases involve using the ESA only. A plot of 
a typical non-repeatable (radiance/gradients) razor effect is 
shown in Figure 1. From this plot, the magnitude of the error 
is seen to be about 200 prad . A Monte Carlo type analysis 
shows that 200-250 prad is about the 3o range fra this error. 
Error angle plots fra the first two cases are shown in Figures 
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2 and 3. With no NR errors in the ESA, the attitude accuracy 
is within 60 prad. With the NR errors in the ESA, this 
accuracy is degraded to about 200 (irad . The large errors in 
the yaw angle estimates are due to filter un-observability. 
The observability of using an ESA combined with gyro 
measurements in a Kalman filter can be shown by using the 
simplifying assumption of a constant coefficient system. The 
state vector in the Kalman filter is given by [1] 


Ax = 


Act 
A b 


(15) 


where Aa is a 3x1 angle error vector (roll, pitch, yaw), and 
A b is a 3x1 gyro-bias error vector. The system error 
equations, state matrix and sensitivity matrix are given by 


Ax = FAx + GAti 
A z = H Ax + Av 


(16a) 


F = 


-[ax] 

®3x3 


~hx3 

^3x3 


(16b) 


— b - - 


o 

0 

1 


(IS) 


which reflects the fact that the spacecraft is Earth -pointing. 
From Equations (16)-(18), the state matrix (F) and 
sensitivity matrix (H) are now constant. 


The observability matrix is given by 


0 = 


H 
HF 
HF 2 


HF 5 


(19) 


which is an 18x6 dimensional matrix. This matrix must be 
rank 6 for the system to be fully observable. However, using 
the system matrices in Equation (16) yields a rank 5 
observability matrix. A singular value decomposition (SVD) 
of the observability matrix can provide an insight to which 
states are observable, as well as the degree of observability. 
The SVD of Equation (19) is given by 


h=[[a( £ )b'x] ; o 3x3 ] (160 

where the angular velocity vector (co) is given by Equation 
(6), and [© x] is the cross product matrix. Therefore, the 
state error angle equations are given by 


USV T = 0 (20) 

where S is an 18x6 diagonal matrix, and U and V are 
unitary matrices with dimensions 18x18 and 6x6, 
respectively. The diagonal elements of the first 6 rows of S 
yield the singular values of the system. These singular values 
yield the degree of observability, which is determined to be 


Act! =— co n Aa 3 — Abj (17a) 

Act 2 = -A b 2 (17b) 

Aa 3 =( 0 „Aa 1 -A ^ (17c) 

The first and third equations show the coupling effects 
between roll and yaw. The nominal body measurements for 
the ESA are given by 


S = 


1 

1 

1 

1 

727x10 

0 


-5 


( 21 ) 


The columns of V shows which states are observable, and 
also show the degree of cross correlated observability in the 
states. This matrix is determined to be 


0 

0 

7.27 xlO -5 
1 
0 
0 


-1 

-2.64 xlO' 5 


0 

0 

0 

0 


-2.64 xlO -5 0 0 0 

1 0 0 0 

0 0 0 1 

0 0 0 -727 xlO' 5 

0 10 0 

0 0-1 0 


( 22 ) 
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The first four columns of V correspond to completely 
observable states. The second and third columns of V 
indicate that the roll and pitch angle states are completely 
observable. Also, there is some correlation between these 
states, shown by the -2.64 x 10~ 5 term. The sixth column of 
V is associated with a singular value of zero. This shows 
that the yaw angle state is not observable. This reflects a 
higher covariance in the yaw angle estimate, as compared to 
the roll and pitch angle covariances (see Figure 4). The 
fourth column of V corresponds to the pitch drift-rate state, 
which is completely observable, since its associated singular 
value is one. The fifth column of V corresponds to the yaw 
drift-rate state, which is weakly observably, since its 

associated singular value is small (i.e., 7.27 xlO” 5 ). 

However, this state is completely decoupled from any other 


state. The first column of V , as well as the sixth column, 
shows the coupling between the yaw angle state and the roll 
drift-rate state (due to quarter-orbit coupling). This indicates 
that the error in this state is attributed to both actual roll rate 
errors and yaw angle errors. Since the yaw angle state is not 
observable, the roll drift-rate errors and yaw angle errors 
cannot be separated. A plot of the gyro drift-rate covariances 
is shown in Figure 5. The error covariance of the roll drift- 
rate state is larger than the yaw drift-rate error covariance. 
This is most likely due to the fact the yaw angle errors cannot 
be separated from the roll drift-rate error. 

Plots of the four remaining cases, which include the 
imager/sounder as an attitude sensor, are shown in Figure 6- 
9. A summary of the results for all cases is shown in Table 1 . 


Table 1 Attitude Errors for Various Sensor Configurations and Error Sources 


Case 

Error Sources 

Roll Errors 

Pitch Errors 

Yaw Errors 

1 

no NR ESA 

60 |irad 

60 (Had 

1x10 s niad 

2 

NR ESA 

200 urad 

200 (Had 

lxlO 5 niad 

3 

no NR ESA, no NR I/S 

60 (Had 

60 urad 

200 firad 

4 

no NR ESA, NR ES 

100(iiad 

100 iirad 

200 (liad 

5 

NR ESA, no NR I/S 

100 urad 

100 urad 

200 niad 

6 

NR ESA, NR I/S 

200 iirad 

200 Miad 

300 (Had 


Since the imager/sounder can measure stars which are off 
nadir, yaw angle information is possible. From Table 1, 
using the imager/sounder as another sensor significantly 
improves the yaw angle estimate. Also, since the magnitude 
of the non-repealable errors is assumed to be approximately 
the same in the ESA and in the imager/sounder assembly, the 
attitude errors are also approximately equal when adding 
these errors to each sensor individually (i.e., case four and 
five). The sixth case involves using both the ESA and 
imager/sounder assembly with non-repeatable errors added to 
each sensor. A purely deterministically found attitude using 
the QUEST method yields errors which are approximately the 
same magnitude as case six (see Figure 10). Therefore, the 
addition of gyros does not seem to significantly improve the 
attitude accuracy. 

In order to possibly estimate the non-repeatable effects in 
the imager/sounder, a colored noise Kalman filter was 
developed. An analysis can be performed by expanding upon 
Farrenkopf s model. The assumed model for the colored- 
noise Farrenkopf analysis is given by [3] 


0 =©£ — T| v —b 


(23a) 


b — x\ u (23b) 

1 = (23c) 

where 0 is the scalar (single-axis) attitude angle, 5> g is the 
gyro output, b is the gyro-drift rate, is the colored-noise 
output, T| v , r\ u , and i] / are zero-mean Gaussian processes 

with standard deviations of a v , a u , and a { , respectively, and 
A[ and B t are the colored-noise system matrices, given by 



r 

o 


(24a) 
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(24b) 


filter. However, an analysis using actual data should be 
performed to investigate the validity of this approach. 



where ©„ is set to orbit rate. Therefore, the full continuous 
system matrix from Equations (23) and (24) is given by 


(25) 

with the state-transition matrix of F denoted by 0 The 
measurement model is given by 


'0 -1 0 O' 

0 0 0 0 

“ 0 0 0 1 

0 0 -co^ 0 



where v is a zero-mean Gaussian process with covariance r , 
and H is given by 

tf = [l o 1 0] (27) 

Equations (26) and (27) show that the colored-noise is added 
to the measurement. The state-noise covariance matrix can 
be computed as 

l/2c 2 Ar 2 0 0 

°** ° ° (28) 

0 0 0 

0 0 c 2 Ar_ 

where Ar is the sampling interval. The steady-state error 
covariance just subsequent to an update is given by 

P = <bP<t> T -PH t \h PH r + rj~ 1 HP + Q (29) 

which can be solved using an eigenvector decomposition of 
the Hamiltonian matrix, where O is the state t ransi tion 
matrix of Equation (25). 

The standard deviation of the colored-noise input varies 
from c i =1x10 2 (jrad to o ; = 1 x 10 ^ prad , which 
corresponds to a colored-noise magnitude ranging from 
15jirad to 1800 prad (these are 3c values). This colored- 
noise output simulates the non-repeatable effect in the 
imager/sounder assembly. A plot of the steady-state colored- 
noise attitude accuracy is shown in Figure 11. Note that the 
standard Farrenkopf analysis with no colored-noise gives an 
attitude accuracy of 56 fir ad (3c value), which is si milar to 
the results shown in Figure 6. This colored-noise analysis 
shows that using an accurate model for the non-repeatable 
errors can reduce the attitude errors when using a Kalman 


c 2 A/ + V3c 2 Ar 
-l/2c 2 Ar 2 
0 
0 


TypicaJ Non-RepeataWe Error for ESA Measurement* 



Pitch 




Errors 
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Star Tracker 

In this section, the simulation results using a star tracker 
with and without gyros are presented. First, the star tracker 
model and parameters are shown. Then, a covariance 
analysis is presented in order to determine the optimal 
orientation of the star trackers. Next, the availability of 
actual stars for the GOES orbit is shown. Results are then 
presented using QUEST [4] to determine the spacecraft 
attitude. An Enhanced QUEST algorithm is also derived 
which filters sensor noise. Finally, simulation results are 
presented using gyros and a Kalman filter. 

All results shown in this section include the dynamics 
and external disturbance in the spacecraft. The GOES Flight 
Software Dynamics Model implements the GOES AOCE 
firmware emulation FORTRAN code from the SS/L into a six 
degree of freedom dynamics model. The initial model was 
developed to examine the replacement of the ESA with gyros, 
and the current capability was developed to compare with 
actual GOES performance using the ESA. A star tracker and 
star tracker/gyro were also added into the simulation. The 
simulation includes rotating solar array inertia effects with 
fully coupled inertia tensor dynamics, magnetic torquers with 
ideal torque response, and gravity gradient and solar pressure 
disturbances. 


The star tracker can sense up to six stars in an 8° x8° 
field of view with a sampling interval of 0.1 seconds. The 
catalog contains stars which can be sensed up to a 6.0 
magnitude. The star tracker measures the tangent of two 
angles, p a and P 2 , resulting in a body vector given by 




1 

^1 + tan 2 p! +tan 2 P 2 


tanPj 

tanp 2 

1 


(30) 


where the z-axis of the star tracker is along the boresight. 
The star tracker ‘"measurements” are obtained by using 


tanj} f =tanp, +v 5( , i = 1, 2 (31) 


where v s is a zero-mean Gaussian process with a 3a value of 
87.2665 prad (18 arc - sec ). 

Each star tracker must be positioned so that sun 
obtrusions can be avoided at all times. For the GOES orbit, 
and available sun shade for the star tracker, the minimum 
exclusion area (allowing for a 3° safety margin) is from 55° 
North to 55° South of the Nadir vector. For the single star 
tracker case, the 55° orientation produces the following order 
for knowledge accuracy: (1) roll angle (i.e., about the 
spacecraft’s x-axis) is known most accurately, then (2) yaw 
angle (i.e., about the spacecraft’s y-axis), and (3) pitch angle 
(i.e., about the spacecraft’s y-axis) being the least accurate. 


The roll is determined to be most accurate since the star 
tracker is perpendicular to this spacecraft’s x-axis. Pitch 
accuracy cannot be improved since the 55° star tracker 
position leads to the y-axis being the least “orthogonal” axis 
with respect to the tracker boresight. 


For the two tracker case, a covariance analysis was 
performed in order to determine the optimal orientation. 
Assuming that each star tracker measures one star for 
simplicity, the error covariance matrix is given by [5] 



+ b 2 bl xb 2 ){b l xb 2 ) T 


(32) 


where 0 is the measurement error standard deviation, and b x 
and b 2 are measurement vectors of each star. For a North- 
South configuration, these measurement vectors are given by 


0 


'O' 

sin55° 

= 

s 

cos 55° 


c 


(33a) 


0 


'O' 

-sin55° 

= 

-5 

_ cos55° 


C 


(33b) 


Using Equation (33), the covariance in Equation (32) 
becomes 


Colored-Noise Kalman Filter Using the Imaget/Sounder Only 



Figure 11 Steady-State Colored-Noise Kalman Filter 
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p=- 


1 

0 


c 

0 0 


0 0 
4 o 


(34) 


The next configuration studied was to place the both star 
trackers 55° North (or South) from Nadir and separated by an 
angle fr , The measurement vectors for this case are given by 



(35) 


where s = sinfr, and c = cos 0 . The covariance matrix in 
Equation (32) for this case is given 


P = 


l(c 4 s 2 c 2 + c 2 s 2 s 2 ) 


2-2 

CS 

0 

0 


s 2 +c 4 s 2 c 2 


esc -c 3 s 2 c s 


csc-c i s 2 cs c 2 c 2 +c 2 s 2 s 2 


(36) 


In order to determine the optimal separation angle, a cost 
function involving roll and pitch errors (i.e., allowing for 
relaxed yaw error conditions) is defined, given by 


m= 




s 2 c 2 + c 2 s 2 s r 


( c 2 s 2 +s 2 + c 4 s 2 c 2 


) (37) 


Minimizing this cost function with respect to ~d leads to the 
optimal separation given by ii= 90° . Therefore, the 
covariance matrix in Equation (36) becomes 


1 


P 


2 



0 


0 

,2 

0 


0 

0 

1 


(38) 


Equation (38) shows that the yaw angle contains the smallest 
error, even though yaw was relaxed for the optimal separation 
angle. Therefore, comparing Equation (34) and Equation 
(38) leads to the conclusion that the optimal location for the 
two tracker case is given by one tracker 55° North and one 
tracker 55° South from Nadir. 


Simulation Results 

Figure 12 shows the actual number of stars within the 
North pointing tracker field of view. There is always a 
minimum of 2 stars, except for the interval from 2.15 to 
2.283 hours. A star with quality 1, but with a magnitude of 
6.256, was added in this interval for the QUEST solution. 
Figure 13 shows the number of stars within the South 
pointing tracker field of view. A quality 2 star (5.137 
magnitude) from the interval 15.45 to 15.483 hours, and 
another quality 1 star (6.138 magnitude) can be added to 
insure a minimum of two stars. This was not done to the 
South tracker catalog, since the North tracker was used for 
simulations involving one tracker. Figure 14 shows the 
combined number of stars for both trackers (without the 
addition of any stars). This shows that a minimum of 4 stars 
is available for the two tracker case. Also, the percentages of 
time in the orbit with the number of stars in the field of view 
are shown by Tables 2 and 3. 


Table 2 North Pointing Star Catalog 


Number of Stars 

Percentage in FOV 

0 

0.0 

1 

0.625 

2 

10.972 

3 

15.625 

4 

27.709 

5 

23.958 

6 

21.111 

Table 3 South Pointing Star Catalog 


Number of Stars 

Percentage in FOV 

0 

0.0 

1 

1.458 

2 

8.056 

3 

20.972 

4 

28.889 

5 

23.272 

6 

17.153 


In this section simulation results using the QUEST and 
Enhanced QUEST algorithms without gyros are presented. 
The QUEST algorithm minimizes the following cost function 
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where A is the attitude matrix, and n is the number of stars 
available. QUEST is a deterministic approach which utilizes 
a point-by-point solution. Therefore, previous measurements 
are not utilized in the attitude solution. This algorithm 
requires at least two star measurements to determine the 
attitude. Therefore, a star is added (as previously described) 
to the single star tracker case. 

In general, the attitude knowledge is determined more 
accurately as the number of star measurements at one time 
increases and/or the separation distance between stars 
increases. This can be seen by the deterministic error 
covariance, given by [4] 




(40) 


*=l 


value which is too large adds too much measurement noise, 
and tends to neglect model corrections. A value of a = 0.05 
was determined to be optimal. The EQA covariance is 
derived by re-writing Equation (41) as 


ijt+i 





(42a) 



(42b) 


where ® denotes quaternion multiplication (see [1]). The 
QUEST determined quaternion is written as 


£*+i 



®£* + i 


(43) 


Figure 15 shows the attitude errors from QUEST determined 
attitude using a single (North) star tracker. Note the large 
errors about 2 hours into the simulation, which corresponds 
to the point where the star availability is mostly only 2 stars. 
Figure 15 also shows the attitude errors along with the 3a 
bounds from Equation (40). This shows excellent agreement 
between theory and simulation. Figure 16 shows the attitude 
errors using both star trackers. Figure 17 also shows the 3a 
bounds for the two star tracker case. This shows the 
significant improvement in attitude knowledge by using two 
trackers. 

In order to further improve the attitude accuracy, an 
Enhanced QUEST algorithm (EQA) was developed. This is a 
simple first-order Kalman filter which combines a propagated 
model with the QUEST determined attitudes. Since gyros are 
not used for this case, the angular velocity is assumed to be 
perfect (i.e., given by Equation (6)). This assumption is not 
true, since external disturbances, and control and sensor 
errors are present in the actual system. Typical control errors 
using the ESA are shown in Figure 17. This shows the large 
errors and dynamic coupling in the roll/yaw axis. The EQA 
is given by 


i k+ , (“) = expj^Q(ffl) Af J q k (+) 

(41a) 

£*(+) = ( 1 -a)£ Jt (-) +a £* 

(41b) 


where At =0.1 seconds, q_ k is the QUEST determined 
attitude at time t k , and a is a scalar gain. This gain was 

determined by minimizing the attitude errors from the 
simulated runs. A value which is too small adds too much 
model corrections, and tends to neglect measurements. A 


where q fc+1 is the true quaternion, and is a three 

component error vector. Substituting Equation (43) into 
Equation (42a), and post-multiplying both sides of the 
resulting equation by yields 


q ~ 

2*+ 1 


'8 q 
Ijt+i 


= O 4 | jt 0£-^(l-a)+a 


8 q 

!*+i 


(44) 


Using a first-order approximation yields the following 
covariance covariance for the EQA 

p m + i = (l-«) 2 ®3 p n k *3 * (45) 

where 0 3 is the state transition matrix of [o> x] . Since this 
matrix is constant and nearly the identity matrix, the diagonal 
elements of Equation (45) approach the following steady-state 

p n k “ T^a P *** (46) 

Figure 18 shows the attitude enrors and bounds from Equation 
(46) using one star tracker and the EQA. Comparing Figure 
18 with Figure 15 shows a significant improvement using the 
EQA. Figure 19 shows the attitude errors using two trackers 
and the EQA. Comparing Figure 19 with Figure 16 again 
shows a significant improvement using the EQA. 

In this section, the results using gyros and a Kalman filter are 
presented. Two gyro cases are simulated. The first case 
involves the utilization of the DRIRU-EL The second case 
involves the utilization of the of the HRG. The parameters 
for both gyros are summarized in Table 4. 
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Table 4 Gyro Parameters 


Parameters 

DRIRU-n 

HRG 

C u (white noise) 
C v (random walk) 

2.15X10 -4 prad/sec 3 / 2 
0.206 jirad/ sec 1 / 2 

L55X10 -4 prad/sec 3 / 2 
1.6 urad/ sec 1 / 2 


The gyro model is shown by Equations (1) and (2). The 
relative performance of the attitude estimation can be found 
by numerically iterating the Kalman filter equations to steady 
state, but Farrenkopf [3] obtained analytic solutions for the 
case when the three attitude error angles are assumed 
decoupled. Farrenkopf s results for the preupdate and 
postupdate attitude error standard deviations, denoted bv 


S v =o v A t 2 /a 


(48d) 


In the limiting case of very frequent updates, the preupdate 
and the postupdate attitude error standard deviations both 
approach the continuous-update limit, given by 




l l 

( 

n 

1 


= At 4 o 2 

a 2 +2c„a v Ar 2 


(47a) 



j 


(49) 


CT ( + )=<*(-)^ (47b) 

where 


l 


r+|s u +(YS H +Sv 2 +js u 2 J 2 

(48a) 

1 

Y=[4 + 5 2 +(l/12)5 2 ]2 

(48b) 

3 
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Using the parameters in Table 4 in Equation (49), it was 
determined that the DRIRU-II steady-state error is 
approximately 2.8 times better (i.e., more accurate) than the 
HRG. This is also shown in the simulations. Figures 20 and 
21 show the attitude errors using the HRG for the one tracker 
and two tracker cases, respectively. Figures 22 and 23 show 
the attitude errors using the DRIRU-II for the one tracker and 
two tracker cases, respectively. Comparing Figure 20 to 
Figure 22, and Figure 21 to Figure 23, it is seen that the 
DRIRU-II is approximately 2 to 3 more accurate for the 
attitude knowledge than using the HRG. Results for the cases 
without gyros and cases with gyros are shown in Table 5 and 
Table 6, respectively. 


Table 5 Attitude Error Results Without ft vrns 


Cases Simulated 

Roll Error 
(|irad) 3a 

Pitch Error 
(|irad) 3a 

Yaw Error 
(jirad) 3a 

QUEST (1 Tracker) 

60 

1250 

900 

QUEST (2 Trackers) 

35 

70 

50 

EQA (1 Tracker) 

12 

225 

175 

EQA (2 Trackers) 

6 

10 

8 


161 
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Conclusions 

This study provided some insightful infonnation for 
using gyros on the GOES spacecraft. It was determined that 
die gyros do not significandy reduce the non-repeatable errors 
in the ESA. This was shown by comparing Figure 9 with 
Figure 10. Since the relative error is approximately equal in 
these two plots, we conclude that the utilization of on-board 
gyros does not significandy improve performance. Also, 
using gyros does not provide any observability in the yaw 
angle estimate, when using the ESA. 

The star tracker simulation results show a significant 
improvement over the ESA attitude knowledge errors. The 
greatest improvements were showing using either: (1) two 
trackers with the EQA, or (2) one tracker and a DRIRU-II 
type gyro, and (3) two trackers and either an HRG type gyro 
or a higher quality gyro such as the DRIRU-n. Adding gyros 
to the spacecraft is the most ideal case since the filter 
bandwidth is larger than the EQA filter bandwidth (i.e., the 
Kalman filter with gyros can sense higher frequency 
spacecraft motions than an EQA). The utilization of on- 
board gyros may also improve the pointing accuracy, since 
the controller bandwidth may be increased. 
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